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Abstract 

A multivariate analysis was performed for populations of Echinocereus (section Triglochidiatus) 
to facilitate the taxonomic circumscription of E. arizonicus subsp. arizonicus. Twenty-one 
morphological characters for 16 populations evidenced the validity of at least two subspecific 
taxa within E. arizonicus'. E. arizonicus subsp. arizonicus and E. arizonicus subsp. nigrihorridispinus. 
Principle components analysis indicated that stem characters were most diagnostic in defining two 
distinct groups of populations, each including the type locality of one of the two subspecies. 
Unweighted pair group method with arithmetic mean (UPGMA) clustered populations of the two 
subspecies apart from one another and from those of the outgroup, E. triglochidiatus subsp. 
mojavensis. For most measured characters, means differed significantly between the two subspecies. 
Discriminant analysis correctly classified 97.0% for individuals of E. arizonicus subsp. arizonicus and 
94.7% for individuals of E. arizonicus subsp. nigrihorridispinus , compared to an overall 97.8% 
correct classification of individuals for all perfect-flowered taxa of section Triglochidiatus 
investigated. 

Key Words: Echinocereus arizonicus , section Triglochidiatus , multivariate analysis, morphological 
characters. 


Taxonomy within section Triglochidiatus H. 
Bravo of Echinocereus Engelm. in F. A. Wislize- 
nus has vacillated dramatically over the past few 
decades. The present study evaluates the circum¬ 
scription of E. arizonicus Rose ex Orcutt subsp. 
arizonicus using multivariate techniques to com¬ 
pare morphological variation within populations 
to that among populations. It is the second 
portion of an ongoing phenetic analysis of section 
Triglochidiatus , the first of which addressed the 
recently discovered species, E. yavapaiensis M. A. 
Baker, and summarized the current knowledge of 
polyploidy within the section (Baker 2006). The 
taxonomic status of E. arizonicus subsp. arizoni¬ 
cus , which is federally listed as Endangered, has 
importance for conservation efforts because of 
threats from mining and other human related 
impacts. 

Until recently, the most widely recognized 
treatment of Echinocereus was by Benson 
(1982), who grouped all of the red-flowered 
populations within the United States under 
a single species, E. triglochidiatus Engelm. in F. 
A. Wislizenus. Since Benson’s treatment, other 
authors reported on the section (Taylor 1985; 
Ferguson 1989) but did not use biosystematic 
approaches. Recent taxonomic, cytological, and 
floral investigations have led specialists to sepa¬ 
rate Benson’s North American E. triglochidiatus 
into at least five species, E. arizonicus , E. 
coccineus Engelm. in F. A. Wislizenus, E. 
santaritensis W. Blum & Rutow, E. triglochidia¬ 


tus , and E. yavapaiensis M. A. Baker (Hoffman 
1992; Blum et al. 1998; Zimmerman & Parfitt 
2003; Baker 2006). The primary rationale for 
splitting E. triglochidiatus into several species is 
the occurrence of polyploidy correlated with 
morphology, geographic distribution, and floral 
dimorphism (Table 1). Echinocereus arizonicus 
represents smooth-spined, diploid, perfect-flow¬ 
ered populations from the Sonoran-Chihuahuan 
Desert interface; E. triglochidiatus represents 
papillate- [E. triglochidiatus subsp. mojavensis 
(Engelm. & Bigelow) W. Blum & Michael Lange] 
or angular-spined (E. triglochidiatus subsp. tri¬ 
glochidiatus ), diploid, perfect-flowered, popula¬ 
tions from the Mojave Desert, California east to 
northern New Mexico and north into Utah and 
Colorado; E. coccineus represents smooth to 
papillate-spined, tetraploid, florally dimorphic 
populations from southern Utah south into the 
mountains of Arizona and from southern Colo¬ 
rado south into Texas and southern Chihuahua; 
E. santaritensis W. Blum & Rutow represents 
tetraploid, per feet-flowered, populations from 
southern Arizona, and E. yavapaiensis hexaploid, 
florally dimorphic, populations from central 
Arizona. Several taxa endemic to Mexico remain 
poorly understood. For a synopsis of chromo¬ 
some numbers in Echinocereus , section Triglochi¬ 
diatus, see Cota and Philbrick (1994) and Baker 
(2006). 

Nomenclature herein follows that of Blum et 
al. (1998). Three taxa are recognized within E. 
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Table 1. Current Taxonomic Status of Echinocereus , Section Triglochidiatus for Populations in 
the United States as Recognized by Zimmerman & Parfitt (2003) and Baker (2006). 


Taxon 

Spine surface 

Ploidy level 

Flowers 

Status in the current 
multivariate analysis 

E. arizonicus 

smooth-spined 

diploid 

perfect 

Analyzed 

E. triglochidiatus subsp. papillate-spined 

diploid 

perfect 

Analyzed as an outgroup 

mojavensis 





E. triglochidiatus subsp. angular-spined 

diploid 

perfect 

Analyzed as an outgroup 

triglochidiatus 





E. santaritensis 

smooth-spined 

tetraploid 

perfect 

Analyzed as an outgroup 

E. coccineus 

smooth to papillate-spined 

tetraploid 

dimorphic 

Not included, see Baker 





2006 

E. yavapaiensis 

smooth-spined 

hexaploid 

dimorphic 

Not included, see Baker 





2006 

arizonicus , E. arizonicus subsp. arizonicus , E. 

Methods 


arizonicus subsp. nigrihorridispinus W. Blum & 
Rutow, and E. arizonicus subsp. matudae (Bravo- 
Hollis) Rutow. According to Blum et al. (1998), 
E. arizonicus subsp. arizonicus occurs in Cochise, 
Gila, Graham, Pima, and Pinal Counties, Ar¬ 
izona, and is characterized as having 7-13 radial 
spines and 1-4 central spines per areole and stems 
with 8-11 ribs; E. arizonicus subsp. nigrihorridis¬ 
pinus occurs in areas southeast of the distribution 
of E. arizonicus subsp. arizonicus , in southern 
Arizona, southwestern New Mexico, and north¬ 
ern Chihuahua and is characterized as having 10- 
14 radial spines and 3-8 central spines per areole 
and stems with 10-13 ribs; E. arizonicus subsp. 
matudae occurs northwestern Chihuahua and is 
characterized as having 7-11 radial spines and 1- 
4 central spines per areole and stems with 6-8 
ribs. 


Twenty one continuous characters (Table 2) 
were measured for at least 30 mature individuals 
from six populations of E. arizonicus , two of E. 
triglochidiatus subsp. triglochidiatus , four of E. 
triglochidiatus subsp. mojavensis , and three of 
E. santaritensis (Fig. 1, Table 3). Because of 
permitting constraints, stem characters only were 
measured for the only known population of E. 
arizonicus subsp. matudae (Bravo-Hollis) Rutow. 
For statistical purposes, the assumption was 
made that all individuals measured within each 
population belonged within a single taxon. 
Although it is possible that individuals of more 
than one taxon occurred within the vicinity of 
any one study population, there were no loca¬ 
tions where this was apparent from the physiog¬ 
nomy of the individuals present. 


Table 2. List of Characters Used in the Analysis. Except for STEML, three measurements were made for 
each character per individual. The youngest fully mature areoles were chosen for stem measurements. 


Character 


Explanation 


STEML 

STEMDIA 

NRIBS 

DBTWARLS 

NCENTRALS 

NRADIALS 

LRADIALS 

LCENTRALS 

THKNCNTR 

FLRL 

FLR WIDTH 

OUTSDIA 

AXIALL 

LTOUPPER 

NTEPALS 

STYLEDIA 

SPINEL 

STYLEL 

NCTRYDIA 

OVARYL 

STAMENL 


Length of longest stem from ground level 

Average stem diameter 

Average number of stem ribs (costae) 

Average distance between stem areoles 
Average number of central spines per stem areole 
Average number of radial spines per stem areole 
Average length of radial spines per stem areole 
Average length of central spines per stem areole 

Greatest thickness of central spine (just above the swollen base) per stem areole 
Flower length from base of pericarp to tip of longest tepal 
Flower width from tepal tip to tepal tip 

Outside diameter of throat of pericarp at constriction just above ovary 
Axial length of pericarpal pith between ovary and base of pericarp 
Length of pericarp from base to uppermost areole 
Number of tepals, excluding those with areoles at their base 
Diameter of style at midpoint. 

Length of longest within uppermost areole of pericarp 

Length of style from its base to bottom of stigma 

Diameter of nectar chamber 

Length of ovary 

Length of stamens collectively 
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Fig. 1 . Locations of study sites for populations of Echinocereus section Triglochidiatus. Numbers next to symbols 
refer to those of populations in Table 1. 


Principle component analysis (PCA) (Systat7, 
SPPS Software Inc. 2000) was used to assess the 
taxonomic values of characters and to assign 
populations of Arizona and New Mexico popula¬ 
tions of E. arizonicus to either E. arizonicus subsp. 
arizonicns , or E. arizonicus subsp. nigrihorridispi- 
nus. This was done by comparing individuals 
within unknown populations to those of the two 
type localities. In order to assess the taxonomic 
value of groups of variables (primarily stem 
characters vs. flower characters), PCA was first 
performed on all characters of populations of all 
three perfect-flowered species, and then performed 
using flower characters and stem characters alone. 
Varimax rotation was used to improve the in- 
terpretability of the scatter diagrams. The cluster¬ 
ing algorithm unweighted pair group method with 
arithmetic mean (UPGMA) was performed with 
NTSYS® 2.1e (Rohlf 2000) to compare phenetic 
distance among populations of E. arizonicus and 
the outgroup, E. triglochidiatus subsp. mojavensis . 
This outgroup was selected because it was the only 
other diploid taxon for which at least three 
populations had been measured. 


MANOVA (SPSS 10) was used to test the 
assumptions of multivariate statistics and to test 
the significance of characters among taxa. Data 
were transformed, as necessary, to meet multi¬ 
variate assumptions. Not all variables met 
homogeneity of variance assumptions after 
transformations. Also, the assumption of homo¬ 
geneity of covariance matrices could not be met 
(Box’s M test, P < 0.001). The Box’s M test, 
however, is generally too strict with the large 
sample sizes generally necessary for multivariate 
applications of ANOVA (Tabachnick 2001). 

Discriminant analysis (DA; SPSS 10) was used to 
test for the correct classification of individuals within 
their respective taxa. Because of the small sample 
size, individuals of E. arizonicus subsp. matudae 
were not included in the analysis. Nine individuals 
were identified as multivariate outliers (Mahala- 
nobis distance-squares from group means with P 
< 0.001) and were deleted from the analysis. A 
permit (TE-844147) for collecting flowers from 
individuals of E. arizonicus was issued by the U. S. 
Fish and Wildlife Service, Southwest Regional 
Office, Albuqerque, New Mexico. 































































Table 3. Locations of Populations (Type Localities in Bold, Chromosome Determinations from Baker [2006]). 
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Results 

Principle Components Analysis (PCA) 

Flower characters dominated the first two 
components of the PCA analysis that included 
all characters (Table 4). When components one 
and two were plotted together (not shown), 
however, there was a poor grouping of individ¬ 
uals within their respective taxa. Component 
three, which also explained a large percent of the 
total variance, had a stronger loading for stem 
characters. When components one and three were 
plotted together, individuals of the three taxa 
began to resolve but there remained a great deal 
of overlap (Fig. 2). For the PCA using only 
flower characters, none of the components 
adequately grouped individuals into species. 

As with the interspecific analysis, stem char¬ 
acters discriminated individuals among taxa of E. 
arizonicus better than flower characters. A 
scatterplot of the first two components of PCA 
defined two groups that maintained the identity 
of individuals within populations (Fig. 3). One 
group was defined by populations one, two, and 
three, and the other group by populations four, 
five, and six. The first group contained the type 
locality for E. arizonicus subsp. arizonicus (pop¬ 
ulation one) and the second for E. arizonicus 
subsp. nigrihorridispinus (population six). The 
highest loadings in the first component were for 
length of both central and radial spines, number 


of central spines, and stem diameter (Table 5). 
The highest loadings in the second component 
were for the distance between areole, thickness of 
central spine, and stem diameter. 

Descriptive statistics, based on the defined 
populations, are presented in Table 6. Those for 
outgroups are presented in Table 7. 

Unweighted Pair Group Method with Arithmetic 
Mean (UPGMA) 

Results from a UPGMA analysis of stem 
characters for Echinocereus arizonicus with E. 
triglochidiatus subsp. mojavensis as an outgroup 
indicated that all populations were placed cor¬ 
rectly within their respective taxa (Fig. 4). As 
expected, the outgroup (populations 11-14) was 
placed with the greatest phenetic distance with 
respect to populations within E. arizonicus. 
Populations of both E. arizonicus subsp. arizoni¬ 
cus (populations 1—3) and E. arizonicus subsp. 
nigrihorridispinus (populations 4-6) were distinct¬ 
ly grouped and the single population of E. 
arizonicus subsp. matudae (population 7) placed 
basally to either of the other two subspecies. 

Multivariate Analysis of Variance (MANOVA) 

Individuals of E. arizonicus subsp. arizonicus , 
as a group, possessed seven characters with 
distinct means: stem length, stem diameter, 


Table 4. Component Loadings for all Characters of PCA with Varimax Rotation of Populations 

OF E. ARIZONICUS (EXCLUDING E. ARIZONICUS SUBSP. MATUDAE ), E. SANTARITENSIS, AND E. TRIGLOCHIDIATUS. 
For translation of character acronyms, see Table 1. 


Component 



1 

2 

3 

4 

5 

STYLE L 

0.9 

0.0 

-0.1 

0.0 

0.0 

FLRL 

0.9 

-0.2 

0.0 

- 0.0 

0.2 

STAMENL 

0.8 

0.1 

0.1 

-0.1 

-0.3 

LTOUPPER 

0.8 

0.2 

-0.1 

-0.2 

0.3 

NRADIALS 

0.0 

0.9 

0.0 

-0.1 

0.0 

NCENTRALS 

0.1 

0.9 

0.0 

-0.1 

-0.1 

DBTWARLS 

0.1 

-0.8 

0.1 

0.3 

-0.1 

NRIBS 

-0.0 

0.8 

0.0 

0.1 

0.3 

OVARYL 

0.4 

-0.5 

0.1 

-0.1 

0.3 

NCTRYDIA 

-0.1 

0.0 

0.9 

-0.1 

-0.1 

OUTSDIA 

0.0 

0.1 

0.8 

0.1 

0.1 

STYLEDIA 

0.0 

-0.0 

0.6 

- 0.0 

0.0 

LCENTRALS 

-0.1 

0.1 

-0.0 

0.9 

-0.1 

LRADIALS 

-0.1 

-0.4 

- 0.0 

0.8 

0.2 

AXIALL 

0.4 

0.1 

-0.2 

-0.2 

0.7 

NTEPALS 

0.0 

0.0 

0.3 

0.3 

0.6 

STEMDIA 

0.1 

-0.1 

0.2 

-0.2 

0.0 

STEML 

-0.1 

0.3 

0.1 

0.1 

0.2 

THKNCNTR 

- 0.0 

-0.1 

0.2 

0.2 

-0.4 

FLR WIDTH 

0.5 

-0.0 

0.3 

0.2 

0.2 

SPINEL 

0.2 

0.5 

0.1 

0.4 

- 0.0 



Percent of total variance explained 




17.0 

17.8 

10.3 

9.8 

7.2 
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Fig. 2. Scatterplot of components 1 and 3 for PCA with varimax rotation using all characters for individuals of E. 
arizonicus , E. santaritensis , and E. triglochidiatus. 
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Fig. 3. Scatter plot of first two components of PCA with varimax rotation including stem characters for 
populations of E. arizonicus , excluding the one population of E. arizonicus ssp. matudae. 
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Table 5. Component Loadings for Stem 
Characters of PCA with Varimax Rotation of 
Populations of E. arizonicus , Excluding that of 

E. ARIZONICUS SUBSP. MATUDAE. 




Component 

1 

2 

3 

NCENTRALS 

0.8 

-0.1 

-0.1 

LCENTRALS 

0.8 

-0.3 

0.4 

NRADIALS 

0.8 

0.0 

0.0 

LRADIALS 

0.7 

-0.3 

0.4 

DRIBS 

0.6 

-0.1 

-0.3 

STEML 

-0.0 

0.8 

-0.2 

STEMDIA 

-0.4 

0.8 

0.2 

DBTWARLS 

0.0 

0.7 

0.5 

THKNCNTR 

-0.0 

0.1 

0.8 


Percent of total variance explained 


30.7 

20.7 

16.6 


distance between areoles, number of central 
spines, number of radial spines, length of 
pericarp, and style diameter (Table 8). Individu¬ 
als of E. arizonicus subsp. nigrihorridispinus , as 
a group, possessed ten characters with distinct 
means: stem length, number of ribs, distance 
between areoles, number of centrals, number of 
radials, length of radials, flower length, pericarp 
spine length, style length, and ovary length 
(Table 8). 

Discriminant Analysis 

Discriminant analysis correctly classified 
97.8% of the original grouped cases that included 
all of the study taxa (Table 9). Individuals of E. 
arizonicus subsp. arizonicus were classified cor¬ 
rectly 97.0%, with 3.0% incorrectly classified as 
E. arizonicus subsp. nigrihorridispinus. Individuals 
of E. arizonicus subsp. nigrihorridispinus were 
correctly classified 94.7%, with 2.1% incorrectly 
classified as E. arizonicus subsp. arizonicus and 
3.2% as E. santaritensis. A single individual (1%) 
of E. santaritensis was misclassified as E. 
arizonicus subsp. nigrihorridispinus and one as 
E. triglochidiatus subsp. mojavensis. Within E. 
triglochidiatus , only a single individual was 
misclassified, which was between the two sub¬ 
species. All between group correlations were 
highly significant (Table 10). The jackknifed 
classification matrix showed a one percent re¬ 
duction in correct classification among each 
group. 

Discussion 

Phenetic analysis presented herein supports the 
recognition of at least two infraspecific taxa 
within E. arizonicus: E. arizonicus subsp. arizoni¬ 
cus and E. arizonicus subsp. nigrihorridispinus. 
Although PCA shows incomplete interspecific 


resolution among the species sampled, there were 
significant differences in means for several 
morphological characters between E. arizonicus 
subsp. arizonicus and E. arizonicus subsp. nigri¬ 
horridispinus. Morphological differences coupled 
with allopatric geographic distribution have long 
been considered basic criteria for the recognition 
of two separate subspecies (Stebbins 1950; 
Lawrence 1951). Most of the known geographical 
distribution for E. arizonicus is represented by E. 
arizonicus subsp. nigrihorridispinus , which was the 
most variable subspecific taxon within the species 
in terms of the diagnostic stem characters. 
Individuals of E. arizonicus subsp. arizonicus , 
which have a limited geographic range, were 
more variable with respect to most flower 
characters. 

Populations of E. arizonicus subsp. arizonicus 
differed significantly from those of E. arizonicus 
subsp. nigrihorridispinus in the means of most 
characters measured (Table 8). However, as 
shown by PCA, flower characters, as a group, 
were not generally diagnostic within the perfect- 
flowered taxa as a whole and no single flower 
character appeared to be useful in separating 
populations of E. arizonicus subsp. arizonicus 
from those of E. arizonicus subsp. nigrihorridis¬ 
pinus. Although stem characters were shown to 
be more diagnostic than flower characters, the 
former tend to be more affected by factors of age 
and environment. Stem diameter, for example, 
should be avoided as a key character because of 
its correlation with available water. Although age 
is generally a critical factor in stem length, stem 
length is mostly determinate within section 
Triglochidiatus because of its cespitose habit. In 
addition, the effects of age were minimized in the 
present study by the measurement of only mature 
individuals. Diagnostic characters that are less 
affected by age and environment are number of 
ribs, number and length of central spines, and 
number and length of radial spines. Although 
even these characters may be affected by etiola¬ 
tion, no individuals occurring in deep shade were 
included in the present study. 

The geographic ranges, as defined herein, of E. 
arizonicus subsp. arizonicus and E. arizonicus 
subsp. nigrihorridispinus differs from that of 
Blum et al. (1998) in that populations of E. 
arizonicus subsp. arizonicus are more restricted. 

Data herein were not sufficient to adequately 
evaluate the taxonomic circumscription of E. 
arizonicus subsp. matudae. Although individuals 
within the single known population possessed the 
fewest ribs, greatest distance between areoles, and 
fewest central and radial spines in comparison to 
those of the other two subspecies, additional 
populations of E. arizonicus subsp. matudae , 
should be sought and measured in order to 
properly address morphological variation 
throughout its range. 







Table 6. Descriptive Statistics of E. arizonicus. 
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Coefficient 

Fig. 4. Phenogram from UPGMA including populations of E. arizonicus and E. triglochidiatus ssp. mojcivensis. 
Matrix composed of mean values for each character. 


Similarly, data were not adequate for the 
assessment of the taxonomic circumscription of 
the two subspecies of E. triglochidiatus , primarily 
because of the lack of data from the type locality 


for E. triglochidiatus ssp. triglochidiatus. Further¬ 
more, UPGMA suggests that E. triglochidiatus 
subsp. mojcivensis may be composed of more than 
one taxonomically definable group, a western 


Table 8. Selected Homogeneous Subsets from MANOVA Significance (Duncans Multiple Range 
Test) where P > 0.02.; Subset A has the Smallest Mean Value and E has the Largest. For all subsets 
containing a single taxon, p = 1.000. 1 = E. arizonicus subsp. arizonicus , 2 = E. arizonicus subsp. nigrihorridispinus , 
3 — E. santaritensis , 4 = E. triglochidiatus subsp. mojavensis , 5 = E. triglochidiatus subsp. triglochidiatus. See 
Tables 5 and 6 for specific means. 


Subset 



A 

B 

C 

D 

E 

STEML 

5 

3 

4 

2 

1 

STEMDIA 

3 

4, 2 

2, 5 

l 


NRIBS 

5 

1, 4 

2 

3 


DBTWARLS 

3 

2 

1 

4 

5 

NCENTRALS 

5 

4 

1 

3 

2 

NRADIALS 

5 

4 

1 

3 

2 

LRADIALS 

1 

3 

2 

5 

4 

LCENTRALS 

1 

5, 3 

4, 2 



THKNCNTR 

3 

4 

2, 5 

5, 1 


FLRL 

2 

4, 3, 1 

5 



OUTSDIA 

3 

5, 4 

2, 1 



LTOUPPER 

4, 5, 2 

1 

3 



STYLEDIA 

3, 4 

5, 2 

1 



SPINEL 

5 

4 

3, 1 

2 


STYLEL 

2 

4, 1 

5 

3 


OVARYL 

2 

3 

1, 4 

5 


STAMENL 

4 

2, 3 

1, 5 
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Table 9. Predicted Group Membership for DA Among Individuals within Populations Preclassified 

AS E. ARIZONICUS SUBSP. ARIZONICUS , E. ARIZONICUS SUBSP. N1GRIHORRIDISPINUS , E. SANTARITENSIS, E. 
TRIGLOCHIDIATUS SUBSP. MOJAVENSIS , AND E. TRIGLOCHIDIATUS SUBSP. MOJAVENSIS. 97.8% of Original grouped 
cases were correctly classified. 



E. arizonicus 




E. triglochidiatus 



subsp. E. arizonicus subsp. 

E. 

E. triglochidiatus 

subsp. 


Taxon 

arizonicus nigrihorridispinus 

santaritensis 

subsp. mojavensis 

triglochidiatus 

Total 



By number of individuals 




E. arizonicus subsp. 

98 

3 

0 

0 

0 

101 

arizonicus 







E. arizonicus subsp. 

2 

90 

3 

0 

0 

95 

nigrihorridispinus 
E. santaritensis 

0 

1 

102 

1 

0 

104 

E. triglochidiatus 




128 

1 

129 

subsp. 

mojavensis 







E. triglochidiatus 

0 

0 

0 

0 

72 

72 

subsp. 

triglochidiatus 









By percent of individuals 




E. arizonicus subsp. 

97.0 

3.0 

.0 

.0 

.0 

100 
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Cereus mojavensis Engelm. & J. M. Bigelow var. 
zuniensis J. M. Bigelow & Engelm, for which the 
type locality is in Canyon Diablo, east of 
Flagstaff, Arizona. 

Evidence from comparative morphology and 
geographic distribution suggests that tetraploid 
populations within section Triglochidiatus , spe¬ 
cifically E. santaritensis , probably arose from E. 
arizonicus subsp. nigrihorridispinus . The two taxa 
are morphologically similar, sympatric over 
much of their ranges, and are both perfect- 
flowered. 


Note that populations display a high degree of 
morphological variability among individuals and, 
consequently, data from several individuals 
should be averaged for identification. Echinocer- 
ens arizonicus subsp. matudae is included in 
a somewhat preliminary sense in that additional 
data may be needed to address morphological 
variation throughout its range. 

A- Stems of mature individuals with mostly 10 
ribs, central spines averaging 4 in number and 
mostly 4 cm long or longer, radial spines 


Table 10. Between Groups F-matrix (df = 21, 701). P < 0.01 for all Values. Wilks’ Lambda = 0.004, 
Approx. F = 70.2 prob < 0.00001. 
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averaging 10 in number. Arizona: Cochise, 
Graham, Greenlee, Pima, and Santa Cruz 
Counties; New Mexico: Hidalgo and Luna 
counties; Mexico: probable occurrence in 

Sonora and Chihuahua. 

. E. arizonieus subsp. nigrihorridispinus 

AA- Stems of mature individuals with mostly 9 or 
fewer ribs, central spines averaging 3 or fewer 
in number and mostly shorter than ca. 3 cm, 
radial spines averaging 9 or fewer in number. . B 
B- Stems of mature individuals with mostly 9 
ribs, central spines averaging 3 in number, 
radial spines averaging 9 in number and 
mostly shorter than 1.3 cm. Arizona: NE 

Pinal and SW Gila Counties. 

. E. arizonieus subsp. arizonieus 

BB- Stems of mature individuals with mostly 7 
ribs, central spines averaging fewer than 2 in 
number, radial spines averaging fewer than 8 
in number and mostly longer than 2 cm. 

Mexico: Chihuahua. 

. E. arizonieus subsp. matudae 
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